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^ ' Abstract: 

m 

We examine the radiative A — > ^yN transition at the real photon point Q 2 = using the frame- 
work of light-cone QCD sum rules. In particular, the sum rules for the transition form factors 
Gm(0) and Rem are determined up to twist 4. The result for G M (0) agrees with experiment 
within 10% accuracy. The agreement for Rem is also reasonable. In addition, we derive new 
light-cone sum rules for the magnetic moments of nucleons, with a complete account of twist-4 
r-| ■ corrections based on a recent reanalysis of photon distribution amplitudes. 
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1 Introduction 



The wish to understand the constituents of atomic nuclei, the nucleons, has been the driving 
force for a great many experiments and theoretical models. In particular the radiative nucleon- 
A-transition has been in the focus of attention since 1979, when it was shown that a deforma- 
tion of the nucleon-A-system can entail non- vanishing electromagnetic (E2) and Coulomb (C2) 
quadrupole amplitudes [1J. This breaks a selection rule laid down previously, which was derived 
in the first non relativistic quark model discussing the E2 amplitude [2J and which allows only 
magnetic dipole amplitudes (Ml) in the 7*iV — > A transition. 

The fact that the measurement of the electromagnetic properties of the transition can pro- 
vide insights in the deviation of the nucleon or the A from spherical symmetry has resulted in 
numerous experiments covering a large range of accessible values for the photon virtuality Q 2 . 
In the whole region up to ~ 4 GeV 2 the ratios E2/M1 and C2/M1 are found to be small and 
negative, especially \E2/M1\ is smaller than 5% jH |H [51 [6l 0. In the case of real photons, 
which is relevant for this work, the Coulomb quadrupole is known to vanish. Thus precision 
measurements are only available for the ratio E2/M1 ll8ll9l [T0l[TT| . 

On the theoretical side various approaches have been suggested. In quark models with un- 
broken SU (6)-spin-flavour symmetry E2/M1 is predicted to be exactly zero, whereas a broken 
SU (6) symmetry yields values ranging from to —2% [fT2l IT3l IT4TI . Other models, among them 
Skyrme models and the large N c limit of QCD, also find the ratio to be small and negative 
lfi~5l[i~6l . Given that the A decays almost entirely into a nucleon and a pion, it is not surprising, 
that chiral bag models tend to agree well with experimental data ifTTl . In recent years studies us- 
ing chiral effective field theory have been quite popular and yield rather precise results lfT8l[i"9l . 
in addition to that lattice calculations also predict E2/M1 to be around —3% GUI . Only re- 
cently a detailed review summarising various theoretical approaches to the nucleon-A-transition 
has been published ([2T]|. 

On the other hand, the attempts to understand the nucleon-A-transition at the microscopic 
level i.e. in terms of the underlying quark-gluon structure have been less successful. In partic- 
ular, the calculation of Ioffe and Smilga [|22l in the framework of QCD sum rules [23J failed to 
produce acceptable results. A possible reason for this is that the background field technique de- 
veloped in [|22l| (see also f24| for an equivalent approach) is only applicable for the case that the 
participating initial state and final state hadrons have equal masses. Technically, this restriction 
arises because the contribution of interest can only be isolated as the double-pole contribution 
in the hadron momentum. This is the case for e.g. the calculation of nucleon magnetic moments 
which was the primary task of If22l 1241 . but it is not a good approximation for the iV — > A7 
radiative transition. 

The problem of calculating the transitions between hadrons of different mass is known for 
a long time and provided the main motivation for the development of an alternative approach 
E51 l26l . now known as light-cone sum rules (LCSRs). In this technique, an infinite series 
of the "induced condensates" (in the language of Il22l ) is resummed in a function that has the 
physical meaning of a photon distribution amplitude and describes the probability amplitude to 
find a quark and an antiquark in the real photon, with given momentum fractions and at small 
transverse separations. The operator product expansion in light-cone sum rules is organised 
in terms of distribution amplitudes (DAs) of increasing twist. The relevant photon distribution 
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amplitudes were introduced in ll25l and recently studied in more detail in [|27l . This technique 
has been used numerously, see e.g. [|28ll29l [301 13TTI for recent applications of photon DAs. 

In this work we calculate the form factors of the 7^ — > A + transition at Q 2 = using the 
light-cone sum rule formalism. In ll32l the transition has been studied at Q 2 = with a LCSR 
approach similar to the one we use in this work. We use, however, a more recent complete set 
of photon DAs up to twist 4, that also include 3-particle-DAs containing an additional gluon, 
which do influence the final sum rules. These were also used in 11331 . where radiative decays 
of decuplet baryons into octet baryons have been considered and in particular the nucleon-A 
transition was also calculated. A comparison with this work is given in Section 3. As already 
stressed in previous works Oil [351, it is important to choose the Lorentz basis in such a way 
that the unwanted contributions due to the non- vanishing overlap of states with spin 1/2 and 
negative parity with the A interpolating field r/ M can be separated from those of spin 3/2 states 
with positive parity. 

For the calculation we will use a technique based on the background field method, that was 
first used in [1251 and Poll to calculate the radiative £ — ► p'j transition and the nucleon magnetic 
moments. In this work we will also give an update on the LCSR results for the magnetic mo- 
ments. Our results for 'jp — > A + can easily be conferred to 772 — > A by exchanging e u <-> e d 
in the final formulae. 

The present analysis is also fuelled by the results of Refs. [[34l [35l . where nucleon-A-tran- 
sition form factors were calculated for virtual photons. In the both calculations that use different 
(local duality and LCSR, using nucleon distribution amplitudes 071 ) techniques, the magnetic 
transition form factor comes out to be below the data for the momentum transfers below 2 GeV 2 , 
and the reason for this discrepancy is not understood. In order to understand the origin of this 
problem it is imperative to have an alternative calculation for the low Q 2 region. Our calculation 
for Q 2 = provides a step in this direction. 

The presentation is organised as follows. In Section 2 we will consider the N — > 'jN 
transition to calculate nucleon magnetic moments. The next Section deals with the nucleon- A- 
transition. We will not give the details of the calculation, but instead focus on the choice of an 
appropriate Lorentz basis. The sum rules for the magnetic dipole form factor Gm(0) and the 
ratio E2/M1 are discussed in Section 4. 

2 Magnetic moments of nucleons 

In this Section we will examine the nucleon magnetic moments. This is a classical problem 
that provides a test ground for many non-perturbative methods. In particular, the calculation 
of nucleon magnetic moments was the main objective behind the generalisation of QCD sum 
rules in background fields [|22l |24| . The results are in good agreement with experiments. In 
Refs. ll36ll38l the magnetic moments were already calculated using LCSRs in conjunction with 
photon DAs. As a new element, our calculation will use, for the first time, the complete set of 
photon DAsl] and also updated non-perturbative parameters. 

'We will not take into account 4-particle-DAs which are not expected to give rise to numerically relevant con- 
tributions 
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Figure 1: A proton with initial momentum p + q emits a photon with momentum q. 



This simple example serves as an illustration how the calculation is carried out and allows 
us to test its accuracy and the dependence of the results on the various parameters. Especially, 
this procedure provides a check for the numerical values of those non-perturbative parameters 
that are still under discussion, e.g. the magnetic susceptibility of the quark condensate x- 



2.1 Definitions 

The transition matrix element 

(P(p,s)\f(0)\P(p + q,s')) = P^(p) 



2m,, 



P (S '\p + q) (1) 



can conveniently be parametrised in terms of Dirac and Pauli form factors, F\ and F 2 . Hereafter 
P^ s \p) is the proton spinor with momentum p and spin s. The magnetic moment of a nucleon 
can then be defined as 



^ = F 1 (0) + F 2 (0). (2) 

This allows us to take only real photons into account. As -F"(0) = and Ff (0) = 1 are just the 
corresponding charges, it is only necessary to determine F 2 (0). 

The process in FigfTJcan be described by the correlation function 



IP (p, q) = i 2 J d 4 x J d A y e^ +i «y (0| T{ Vp (x) f (y) rj p (0)} |0> . (3) 

Here 

f = ejdj^d + CuU-i^u (4) 

is the electromagnetic current, with e d = —1/3 and e u = 2/3 being the quark charges, is 
the four-polarisation vector of the emitted photon, q ■ = 0. The current 

V (x) = (« a (x)C 7 V(x)) l5lx d c (x)e abc (5) 
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is the usual Ioffe current ||39l for the proton. C is the charge conjugation matrix, a, b, c are colour 
indices and e abc is the three dimensional Levi-Civita symbol. Note that the current has fixed 
isospin 1 /2. The isospin relation between proton and neutron then assures that the formulae for 
the neutron magnetic moments can be obtained from those of the proton by exchanging e u <-> e^. 
Therefore we will only consider the proton. The coupling Ap of the Ioffe current Eq.© to the 
proton is defined by 



<0|t/(0) \P(p,s)) 



A; 



(2tt)' 



r^ (s) (p) • 



By introducing an electromagnetic background field of a plane wave 



\J,V 



J rii/C- 



(6) 



(7) 



it is possible to absorb the emitted photon into the background field. We consider the following 
object 



Hp (p, q) ef 



i d'xe^(0\T{ Vp (x)v p (0)}\0) 1 



(8) 



Here the subscript F indicates that the v.e.v. has to be evaluated in the background field F^ v . 
Expanding the correlator in Eq.© in powers of the background field and taking only the terms 
linear in F pu corresponds to the single photon emission, also described by Eq.©. We will refrain 
from a more detailed presentation of this expansion, but instead would like to refer to HOI for a 
review on the background field method and to Chapter 2 in Ref. 11271 . which is dedicated to the 
expansion of correlation functions in an electromagnetic background field. 



2.2 Expansion of the correlator 

Applying the usual strategy of QCD sum rules, we have to calculate the correlation function 
Eq.® in two different regimes. If the in- and out-going particle states are close to proton mass- 
shell, i.e. p 2 « m 2 P and (p + q) 2 ~ m 2 P , the hadronic representation of Eq.® will be dominated 
by the process p — > pj. 

The contribution T P of the process p — > p7 is then given by 



Tp(p,q) = e 



(A) 



<0| 7^(0) \P(p, s)) (P{p, s) | ^(0) \P{p + q, s)) (P{p + q,s)\ r}(0) |0) 



(m 2 p — p 2 ) (m 2 p — (p + q) 2 ) 
Using Eqs.([T]),© and the spin summation formula for Dirac spinors 

Y j P {s \p)P {s \p)=f + mp 



(9) 



(10) 



this can be written as 



|A P | 2 f + m p 
(2n) i m 2 -p 2 



YPiiQ' 



2m. 



-cr^q^Q 2 



'+0 + m P c (X) 

m 2 v -p\ M 



(11) 
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with pi = p and p 2 = p + q. 
The Lorentz structure 

fo^fq v ef 

is free of contributions due to -Fi(O), already satisfies the Ward identity, as it is proportional to 
f^, and has the highest possible number of momenta p. Therefore it seems advisable to focus 
on structures containing an even number of 7-matrices. 

After reducing to the Dirac basis we will only consider the structure p a p l3 a 13 ^ f a ^, following 
1551 . Then one gets for Eq.(HB 

The dots represent terms of different Lorentz structure. 

In the Euclidean region, where p\ <C and p\ <C 0, one can express the correlation function 
Eq.© in terms of photon distribution amplitudes of increasing twist. To this end we insert the 
expressions for the current 77, Eq.©, into the correlator ©. 



n M1 >, q)ei x) = % Jd 4 x e ipx+iqy (0| T { (u a (x)C-f X u b (x)) ^ x d c (x)- 

■ d°' (0) 7 V 75 (V (0) lvCu h ' (0)) } |0) F e abc e a ' b ' c ' (13) 

Using Wick's theorem the calculation is straightforward. One has, however, to pay attention to 
the fact that we are working with massless quarks in a simultaneous electromagnetic and gluonic 
background field. The electromagnetic field F M „ is just a plane wave, whereas the gluonic field 
G>„ = G^ v t A due to the presence of gluons in the hadron is unknown. The quark propagator 
then adopts the following forml pTfl : 

# r nl W_ 

2tt 2 x 4 L ' 1 16tt 2 x 2 



q(x)q(0) = - 2 4 [x, 0] - 2 2 / du [x, ux] {ufia a/3 + ua a ^} G a(3 {ux) [ux, 0] 



16tt 2 x 2 

here we used the abbreviation 



/ du [x, ux] {ufoap + ucTapft} F a/3 (ux) [ux, 0] + . . . , (14) 
Jo 



[x,y] 



Pexp jz J dt{x- y)^ [e q A>*(tx - ty) + gB»(tx - ty)] J (15) 



for the path-ordered exponent (A^ is the electromagnetic and the gluon field) and a = 1 — a, 
Va G [0,1]. The dots represent terms that will give rise to contributions of twist 5 or higher. As 
we will only consider terms up to twist 4, these are not relevant here. Hence, there are only four 
diagrams that have to be taken into account, see Fig|2l 

It turns out that the diagram in Fig. 2a does not give rise to contributions with an even number 
of 7-matrices and can be neglected. 



2 Note that our sign convention for the electric charge, e = ^/47ra em , follows ll27l and thus differs from 
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diagram a 



diagram b 




diagram c diagram d 



Figure 2: Diagrams up to twist 4. The wiggled and the curled lines represent the coupling to 
the electromagnetic and gluonic background fields. The crosses denote interactions with the 
vacuum. 



After using Fierz identity to decompose (0| <7i(a?)?j(0) \0) F , one can insert the expression for 
the photon DAs (see Appendix |A]) and perform the Fourier transformation. 
We get for diagram b: 



T b P (p,q) 



e<i (qq) 

67T 2 

£d (qq) 



du tp(u) In ( _ » x 

1 upf — up^ 



+ 



<,7T Z 



\ A(u) + B(u 
du - 



o 



—upf 



upl 



p>^ k U + 



(16) 



here the dots represent terms that do not contribute to the structure p a ppa^ f a ^ or that are just 
polynomials in p\ and p\. These will vanish after a subsequent double Borel transformation. 
The calculation of the remaining two diagrams is analogous and yields 



T c P (p,q) 



T d P (p,q) 



e u (qq) 
2tt 2 

ed (qq) 

4.7T 2 

| e d (qq) 

471-2 



du / Da 



Sj(a) 



-a u p\ - a u p\ 



du / T>a- 



L - 2u)T2(a) 
-a u p\ - a u p\ 



p a Pf3^ ( f a( + 



du T>a 



S(a) + 5(a) 
-a u p\ - a u p 2 2 



du Va(l-2u) 



T 2 (a)-2T 3 (a)+T 4 (a) 
-a u p\ - a u p\ 



(17) 
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x p a Pf3 a K f a( + ... , (18) 

where J Da = da q da-g da g 5(1 — a q — a-q — a g ). The functions ip, which is of twist 
2 and T- t , S, S, «S 7 , Tj, A and B, which have twist 4, are defined in AppendijfAl 

2.3 Borel transformation and continuum subtraction 

The sum rule for F 2 (0) can readily be obtained be equating the hadronic result, Eq. (fT2]) . and 
the light-cone expansion, Eqs.(TT6l [T71 [T8l) . As usual, a Borel transformation and a subsequent 
continuum subtraction are necessary to suppress the effects of excited states and extract the 
p — > £>7 ratio. The two independent momenta pi and p 2 allow a double Borel transformation, 
that can be performed using the general formulae 

r(q) ) f r 

-up 2 - M p 2 ) Q J V M 2 + M| 

and 



r2 J : I - p-mi/Ml-mZ/Mj (2 r)\ 

BM ^\(ml-pl)(ml-pt)j ■ (20) 

9 9 Af 2 .M 2 

Here M/ is the Borel parameter corresponding to pf and t = M i +] ^2 • The continuum subtraction 
can be accomplished by a simple set of substitution rules P42l : 



t 3 — t 3 (l-e-^(l + ^ + ^(5oAfjJ (2D 
t 2 — ^fl- e -Wl + ^ (22) 

t — ► t (1 - e~ 5o/t ) (23) 

with 5 being the continuum threshold. 

For the process p — ► ^7 the natural choices are 

M 2 = Ml (24) 

and the threshold Sq coincides with the normal continuum threshold sp for the proton. 
Putting everything together the final sum rule for F 2 (0) takes the following form: 



^(0) 

\A 

ed 



7T 2 m p (qq) p t 



ApI 2 



y^(l/2)t 2 (l - e~ s °/* (l + ^ 



+ 



-— (A(l/2) + B(1/2))J t (1 - e- 5 °/*) 



ed 



1/2 ,1/2 1 



— — / da„ / da 



2 ./o JO 1 — a q - a Q 
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s + s + ^s, 

e d 



(at q , a-g, 1- oc q - otq)t (l - e So/t ) 



1/2 



X 



T 2 - 2T 3 + T 4 - — Tl 

en 



(a q , aq,l-a q - otq)t (l - e 5o/t ) 



(25) 



The first term in Eq.(l25l) gives the leading twist-2 contribution, which was first obtained in 
Ref. ll36ll . The remaining terms are new. 



2.4 Numerical Results 

The asymptotic expression for the photon wave function (p(u) is given by x(/i)6u(l — u) at a 
renormalisation scale /j 2 = 1 GeV 2 , where x(aO is the so called magnetic susceptibility of the 
quark condenstate. It has been argued that the full DA does not differ much from the asymptotic 
expression [|36ll as sum rule calculations showed a small coefficient for the next-to-leading order 
term. Thus, henceforth, we will use ip (1/2) = 3/2x(^)- 

The value of \ is n °t very well known. The first detailed study using QCD sum rules 
yielded x(l/2) = 4.4 GeV" 2 03 El, whereas a local duality approach J2l found x(l/2) « 
3.3 GeV" 2 . The latest estimate for x(n = lGeV) gives a value of 3.15 GeV" 2 G71 and, assum- 
ing asymptotic DAs, leads to <p(l/2) « 4.73GeV" 2 . 

The higher-twist photon DAs are known to next-to-leading order in conformal spin. The 
corresponding expressions are collected in Appendix lAl 

Furthermore, we need the numerical values for the coupling constant |Ap| 2 , the continuum 
threshold sp as well as the Borel window for t. The continuum threshold can be determined 
from the most fundamental sum rule for the nucleon, namely those for the coupling constant. 
We obtain sp ~ 2.25 GeV 2 , see ll46ll . It is advantageous to use the whole sum rule ll47l instead 
of a fixed number for |Ap| 2 . This will decrease the dependence of our LCSR on the value of the 
quark condensate and thus improve stability and reduce errors. We are left with the choice for 
the Borel window. 

The Borel window is determined by two competing requirements. On the one hand, t must 
be large enough to ensure that severing the twist expansion after twist 4 is valid, as the contri- 
butions of twist 5 and 6 are suppressed by an additional factor \ compared to the twist 3 and 4 
contributions. On the other hand, a small Borel parameter is necessary to assure an adequate ex- 
ponential suppression of the continuum and guarantees the validity of the quark-hadron duality. 
This suggests the interval 

lGeV 2 < t < 2 GeV 2 . (26) 

In Figj3]we plotted the sum rule (|25T) for different values of ip(l/2). The comparison with 
the experimental value J45l \i P w 2.793 favours a value y?(l/2) = 5.25 ± 0.15 GeV" 2 . This 
corresponds to = 1 GeV) = 3.5±0.1 GeV -2 , assuming one uses the asymptotic expression 
for the DAs. This agrees rather well with the result from the vector dominance model and the 
latest QCD sum rule result. Note that a larger value of x can still be realized, if the full DA has 
a rather flat shape, this is the case e.g. in the instanton-model of the QCD vacuum [48ll49l. 



8 



3,2 



3,1 



(p(l/2)=6.6 GeV 

<p(l/2)=5.0 GeV" 2 

(p(l/2)=4.73 GeV" : 
exp. value 




t in GeV 



Figure 3: Magnetic moment of the proton from sum rule (|25T) for different values of (p(l/2). 
The hatched line represents the experimental value P31 . 

The magnetic moment of the neutron can, as already stated, be obtained from Eq. ((25T) by 
exchanging e u <-» e d . We plotted the results in Fig. HI As in 0611 . the sum rule prediction for fi N 
is somewhat below the experimental value H51 [in ~ —1.913. 

The experimental result for the neutron magnetic moment //jy can be reproduced by using a 
40% larger value of 0(1/2) that can be achieved if the magnetic susceptibility \ is increased by 
the corresponding amount, or if the photon distribution amplitude is more peaked in the middle 
point as might be suggested by the model for the distribution amplitude for transversely polar- 
ized mesons in ll50l . The stability of the neutron sum rule becomes, however, somewhat worse 
in this case. Also also the agreement of the sum rule prediction for the proton magnetic moment 
is spoiled. Also /xjv turns out to be more sensitive than fj, p to those higher twist corrections that 
are only known to circa 50% accuracy so that the overall error is larger. Worse agreement for 

compared to /ip is therefore no surprise. 

With the standard choice (p(u) ~ (p asy (u) and the latest value = 1 GeV 2 ) = 3.15 ± 
0.3 GeV" 2 ||27l3 we obtain the following results for the magnetic moments 

/jp = 1 + (0.96 ± 0.1)^-2 + (0.72 ± 0.18) to _ 4 

= 2.68 ± 0.28 (27) 

H N =(-1.93 ± 0.2) to _ 2 + (0.88 ± 0.27) to _ 4 

3 X only appears in the combination (qq) x, which has a very weak scale dependence, see e.g. |27|. Thus the 
evolution effect is negligible. 
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-0,4 



-1,6 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



■ ■ ■ cp(l/2)= 6.6 GeV 

- ■ cp(l/2)= 5.0 GeV 2 

(p( 1/2)= 4.73 GeV" : 

exp. value 




j i i i I i i i i I i i i i I i i i i I i i i i_ 

1 1,2 1,4 1,6 1,8 2 

t in GeV 2 



Figure 4: Magnetic moment of the neutron from sum rule (l25l) for different values of (/?(l/2). 
The hatched line represents the experimental value P3Tl . 

= 1.05 ± 0.47. (28) 

Here we also included the different contributions of twist 2 and twist 4 to the magnetic moments 
and how the overall error is distributed. 

In the following we will stick to the above standard choice for (p(u). 



3 The Nucleon-Delta-transition 

We will now expand the technique presented in the previous Section to the case of the nucleon- 
A-transition. We will not go into the details of the calculation, except for the construction 
of a suitable Lorentz basis. This enables us to remove unwanted contributions of transitions 
including final states with isospin 3/2 and spin 1/2. 

It is not necessary to treat the the wy — > A transition separately from the — > A + tran- 
sition, as the final formulae will only differ by the exchange e u <-» e d . Hence, we will only 
consider the proton transition. 

3.1 Definitions 

In order to study the pj — > A + transition using LCSRs, it is convenient to consider the correla- 
tion function corresponding to the diagram in FigJ5J 
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1/(0) 




rf{x) 



Figure 5: Graph corresponding to Eq.(l29l) in coordinate space. 



W u (p,q) =i 2 jd 4 x jd A y e ipx+iqy (Q\T{ri»(x)f(y)f](Q)}\Q) . (29) 

For the current 7? M that creates states with the quantum numbers of the A + , we will follow a 
suggestion by Ioffe [39] and use 

Tf( x ) = [(u a (x)C^u b (x)) d c {x) + 2 (u a (x)C^d b (x)) u c {x)] e abc . (30) 

Note that the current rj, Eq.©, has isospin 1/2, whereas rf has isospin 3/2, therefore only 
the isovector part of the electromagnetic current 



$ 1} = g ( e « ~ ed ) ( u ^v u ~ d ^ d ) 

can induce the ^7 — > A + transition ||34ll . 

We will also use a correlation function in an electromagnetic background field F^ u 

n%(p,q)eW=i fd t xe i *°(0\T{7f(x)fj(0)}\0)F ■ 



(31) 



(32) 



The quantities, which can be measured experimentally, are the magnetic dipole transition 
form factor Gm(0) and the electric quadrupole transition form factor Ge(0)- They are given by 

ED 



G M (0) 



mp 



3 (m P + m A ) 



(3mA + mp) (r/iA + mp 



Gi(0) 
mA 



+ (m A - mp) G 2 (0) 



G E (0) 



3 (m P + 771a) 



m A — m, 



Gi(0) 
m A 



G 2 (0) 



and 



EA4 



gg(0) 
" G M (0) 

where Gi(0),i = 1, 2 are form factors defined in Eq.(l39l). 



(33) 
(34) 

(35) 
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3.2 Contribution of jry — > A + to the correlation function 

The contribution of the p'j — ► A + transition to the correlator (l29l) is given by 

1 1 



A " m 2 ^ — p 2 mp — (p + g) 2 ^—j 



53<0|^(0)|A(p >a )> 



(A(p, s) | f(0) \P(p + q, s')) (P(p + q, s>) 1 7/(0) |0) e 



(A) 



The matrix element 



^(0)|A(p,s)) 



Aa iA«(p) 



(2*0 



(36) 



(37) 



is parametrised via the coupling constant Aa- Here A^\p) is a Rarita-Schwinger spinor for the 
A. Inserting the general expression for the transition matrix element 



(A"(p, s) | j„(0) |P(p + q, s')) = A M (p, s)r^ T5 P(p + g, s') , 
where the vertex defines the three form factors Gi(Q 2 ), G^Q 2 ) and Gs(Q 2 ), 

I> = Gi (Q 2 ) (g^ - q^ u ) + G 2 (Q 2 ) (g^q ■ [p + ~J - g M 

+ C 3 (Q 2 ) (g M g^ - g 2 ^) , 
in Eq.(l36l) and using the spin summation formula 

E A ? ) (p)4 s) w = -( 



(38) 



m A 



^ 3 W 3m 2 



+ 



3mA 



(39) 



(40) 



the contribution T^ey can be written as 



A " ' (m 2 A -p 2 )(m 2 P -(p + q) 2 ) 



- m A )x 



3 3mi 



+ 



3m,A 



(41) 



However, it is known that the correlation function is plagued by transitions pj — > iV* 11341 
where \N*) is a isospin-| spin-| state with negative parity. It is advantageous to use the Lorentz 
structures that do not receive contributions due to transitions to spin-| states. To find these we 
consider the overlap of \N*} with jf, which is defined 11351 as 



% (0)\N*(p,s)) 



(27T) 



(42) 



where \ N * is the coupling and the mass of the spin- 1/2 state. The spinor N*^ s \p) satisfies 
the Dirac equation (jf> — m*)iV*w(p) = and equation (fTOl) . Using the general decomposition 
of the electromagnetic transition matrix element 



(N*(p, s)\f(0) \N(p + q, s')) =N* (s \p) [( 7 Y - <tf) P, NN \Q 
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-ia™q a F 2 NN *(Q 2 )] l5 N^ (p + q), (43) 
we can write the unwanted contribution to IP" as 

T ( r /2) =r~ 2 2 ^2 N < 7 ,> >7^ - V)(|> - m^)* 

(m 2 N , -p 2 )(m z N - {p + q) 2 ) 

[{l v q 2 - fq v ) F? N *(Q 2 ) - m^F^Q 2 )] 75 (f> + f - m N ) . (44) 

In Refs. ll3"4ll35l it has been shown that it is possible to disentangle the contribution to IP" due 
to Eq. ([44|) and Eq. (|4Tj) by a specific choice for the Lorentz basis. As our kinematics are different 
from those in Ref. ll35l . we cannot use the same basis, however. 



3.3 Lorentz basis 

The correlator Eq.(|29l) satisfies two independent constraints, which have to be taken into account 
when constructing a suitable Lorentz basis: 

• the transversality condition q^IP" = 

• the Rarita-Schwinger condition 7 Ai rP l/ = . 

As the transversality condition is automatically fulfilled by Lorentz structures proportional to 
F^ v and its derivatives, it is convenient to construct a basis from these structures. This yields 20 
different Lorentz structures and the QED Bianchi Identity eliminates four thereof. The Rarita- 
Schwinger condition provides four additional constraints reducing the number of independent 
Lorentz structures to a mere 12: 



Tit 


= qp Yifii 


5 - ep 7 M ^7 5 + 4 (qp jffa - ep p^7 5 ) 


n 2 


= p 2 (qp 7 M ; 


^75 - ep 7^75) + 4 (qp p M ^7s - ep puffin) 


n 3 


= qp 7^75 


- ep 7^75 + 2 (p"#7s) - - (Y^ 5 ) 


n 4 


= 4 (jffyf) 


'5) - P 2 (7^75) + 2 (qp 7^75 - ep 7^75) 




= qp (i^it, 


/ 5 ) - 4 (qp e M ^7 5 - ep 9^75) 




= 7^75 - 


2 (e^ 7 5 - ?77 5 ) 


n 7 


= qp 7^75 


- ep 7^75 - 4 (qp e M 7 5 - ep g M 7 5 ) 


n 8 


= 4 (qp 


^75 - ep q^pys) - qp (7^7*75) 


n 9 


= 2 (q^ffs 


- e M ^7s) + Y^fls 




= 4 (qp e^f 


75 - ep ^75) - qp 7^75 + e P 7 M ^75 


n u 


= g M ^75 











Here we inserted the explicit expression for f^ v , see Eq.©, which is advantageous for the further 
calculation. 
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3.3.1 The p7 — > A + contribution 



The contribution of the nucleon-A-transition to the correlation function, Eq.(l29l). can be ob- 
tained readily from Eq.(l36l). Using the basis (1451) . the result has the form: 



T£ u (p, q)e 



(A) 



- AaA 



1 



X 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



(m 2 A — p 2 ) (m 2 p — (p + qf 



(p + g) 2 Gi(0) m P G 2 (0) (p + q) 2 G 2 (0) 
96m A vr 4 192tt 4 192m A 7r 4 

mpGi(O) G 2 (0) m P G 2 (0) 



96m A vr 4 



+ 



192tt 4 192m A vr 4 
mpGi(O) p 2 G 1 (0) , 2gpG 2 (0) 2^(0) 



+ 



48tt 4 48m A vr 4 192vr 4 48m A 7r 4 
G?i(0) m P G? 1 (0) 2^(0) 2gj9G 2 (0) 



96tt 4 96m A vr 4 192m A vr 4 384m A vr 4 
Gi(0) , m A G 2 (0) 



■IIa 



+ 



32vr 4 64vr 4 
m P m A Gi(0) p 2 Gi(0) 2gpG 1 (0) 



32tT 4 



32tt 4 



32tt 4 



mpGi(O) mpm A G 2 (0) p 2 G 2 (0) 2qpG 2 (0) 



32vr 4 
G 2 (0) 



48vr 4 



48vr 4 



48tt 4 



64vr 4 

/mpGi(O) m A Gi(0) 



+ 



32tt 4 32tt 4 
Gi(0) mpG 2 (0) m A G 2 (0) 



32tt 4 

gi(0) 
16tt 4 



48vr 4 



48vr 4 



■ n 
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3.3.2 The iV7 -> N* contribution 

Analogously one finds for the contribution of the J p = | states to Eq.(|29]): 



T^ 2 {p,q)e ( u 



(A) 



Xn*Xn- 



m 1 — p 2 ) ( m p — (p + q) 2 ) 

4m P F 2 *(0) ■ - 4F 2 *(0) • TZ 2 
f (4p 2 F 2 *(0) - 4m P m»F*(0)) • ft 3 



Tin 



(46) 
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+ (2m*F 2 *(0) - 2m P F 2 *(0)) • TZ 4 



(47) 



We see that only the Lorentz structures TZ\, IZ2, and 7^ 4 receive spin-l/2-contributions. In 
addition to that, the coefficient of the linear combination 

Tlx + m P n 2 (48) 

vanishes for the A/7 — > N* transition. Thus, we can choose among 9 structures which are 
suitable for the two sum rules for Gi(0) and G 2 (0). 

3.4 Sum rules and numerical results 

The light-cone expansion of the correlation function is completely analogous to Section [2l even 
the Feynman diagrams are identical, see Figj2l The calculation can be simplified further if the 
isospin relation requiring an overall factor (e u — e^) is taken into account. It is then sufficient 
to calculate only terms proportional to e d , which are simpler than those proportional to e u . The 
lengthy results for the various diagrams can be found in Appendix [Bj 

Before writing down sum rules for Gi(0) and G 2 (0), it is advisable to identify those Lorentz 
structures promising the most reliable results. We will use two criteria: 

1. the structure is free of spin- 1/2 contributions 

2. the structure has the highest possible power of the momentum p. 
These demands are fulfilled by 

^5 = QP (7^75) - 4 (qp e^7 5 - ep g^7 5 ) , 
TZ 8 = A (qp e^7 5 - ep 5^75) - qp (t m ^|*7b) » 

and 

K u = q^fls ■ 

Upon collecting the corresponding terms from Eqs. (|4^1) ,(|92l,(l93T),(l94l) and (1951 ), one can easily 
assemble the sum rules corresponding to the three structures. 

The necessary Borel transformation can again be performed with Eqs.(fT9l). (l20l) . There is, 
however, a subtlety: the choice of the ratio of the two Borel parameters Mf and Mf . If we 
were able to calculate the correlation functions exactly, the dependence on the Borel parameters 
would vanish, as they are not physical quantities. Our calculation is, however, approximate, 
and the approximation is rather crude so that the difference in the mass scales in the A and the 
nucleon channels is not reproduced by the sum rules. This can be checked by the calculation 
of the masses in the two momentum channels using standard techniques. This mass difference 
between the nucleon and the A is large, of order 300 MeV, and it has to be taken into account. 
Neglecting this difference, from our point of view, is the main reason why the Ioffe-Smilga sum 
rule for the nucleon A magnetic transition did not produce acceptable results. In the approach 
that we are using there is a possibility to take into account the mass difference because the two 
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momenta pi and pi alias Borel parameters M\ and Mf are taken as independent variables, and 
the sum rule can be "repaired" by taking the ratio of Borel parameters at a fixed value 



Ml 



m. 



m 



(49) 



N 



Since this ratio determines the momentum fraction ratio, at which e.g. the photon wave function 
ip(u) is evaluated, this choice shifts u away from the centre 1/2. We will follow this strategy, 
which was advocated in lf25l for the calculation of the asymmetry in the £ — ► jry decay. A 
similar trick is often used in the calculation of SU(3) flavour symmetry breaking effects in the 
sum rule method. 

The final sum rules for Gi(0) and ^2(0) are given by 



Gi(0) = 16tt s 



ed) (qq) 



,M 2 /t 



-So/t 



1 + 



So 



(50) 



and 



-S /t 



So 

t 



x 9 



(51) 



and 



Gi(0) + ^G 2 (0) 



l7T 2 {e u -e d ) {qq) cM 2 /t 



X A X } 



vB(v) 



t ( 1 _ e -SoM +I 



(52) 



with 



M 



2 2m\rrip 
ml + m\ 



So 



2sls 



2 



5 A + S p 

r'2 



Mf 



Mi + M{ 



(53) 
(54) 
(55) 



Xi and X 3 correspond to lengthy contributions of twist 4, whereas X 2 is of twist 3. The full 
expressions for X 1; X 2 and X 3 can be found in Appendix O The magnetic dipole form factor 
Gm(0) and the electric quadrupole form factor Ge(0) can be obtained via Eqs. (|33l) . d34l) . 

In the following analysis we use the asymptotic expression for the leading photon wave 
function (p(u), which yielded a reasonable result for the proton magnetic moment in Section 2. 
The continuum threshold for the A, sa ~ 3.0 GeV 2 , can be determined from the sum rule for 
the coupling |A A | [46]. By using the whole sum rule expressions for A a and Xp, the stability of 
the sum rule can again be improved. The same Borel window as in Eq.(T26l) is used. 

In Figj6]our result for Ga/(0) is shown. In addition to our sum rule, we also plotted the 
results for a slightly changed ratio of the Borel parameters (±15%). The result implies that the 
sum rules for Gm(0) are rather stable with respect to a variation of v, if the Borel parameters 
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— v=0.633 






■ ■ ■ -15% 






— ■ +15% 






w- exp. value 



























■ ■ ■ -15% 






— v=0.633 






- ■ +15% 






w/. exp. value 











Figure 6: The magnetic form factor Gm(0) of the pj — ► A + transition. The solid black curve 
shows the result of the corresponding sum rule, the dashed and dotted curves arise from a varia- 
tion of the ratio of Borel parameters by 15%. The experimental value [|52l including uncertainties 
is given by the hatched region. The left panel shows the result using G 2 (0) determined from sum 
rule Eq.dFTT). the right panel using G 2 (0) determined from Eq.(l52l). 



are chosen as in (|49b . The numerical values obtained if one determines G 2 (0) from Eq.dFTT) or 
Eq.(l52l). respectively, 

G M (0) = (1-55 ± 0.15) to _ 2 + (-0.04 ± 0.01) to _ 3 + (1.08 ± 0.2) to _ 4 

= 2.59 ± 0.36 (56) 

G M (0) = (1.55 ± 0.15) tw _ 2 + (1.15 ± 0.12)^-4 

= 2.70 ± 0.27 (57) 

are close together and agree rather well with experiment ll52l 

G M (0) = 3.02 ±0.03 . (58) 

Note that the error in the leading twist contribution stems almost exclusively from the uncertainty 
of the value of the magnetic susceptibility x- 

Our estimates for the ratio Rem that can be obtained by using Eq.(T35l) are shown in FigJTl 
The results are 

Rem(0) = (-7.6% ± 0.1%) tw _ 2 + (1-15% ± 0.7%) to „ 3 + (0.05% ± 0.02%) to _ 4 

w -(6.4 ±0.8)% (59) 

Rem(0) = (-6.3% ± 0.35%) to _ 2 + (3.7% ± 0.25%) to _ 4 

w -(2.8 ±0.6)%, (60) 

where G^O) is determined from the sum rule Eq.dFTI) and Eq.(l52l). respectively. These values 
have to be compared to 

Rem(0) = (-2.5 ± 0.5)% , (61) 
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Figure 7: The ratio Rem- The identification of the panels follows Figj6l The experimental value 
is from R5Tl . The left panel shows the result using G 2 (0) determined from sum rule Eq.dFTT). the 
right panel using G 2 (0) determined from Eq.(l52l). 



given by the Particle Date Group ll45ll . Although the result (|6Q|) is closer to experimental data, 
it is less reliable, as Eq.(l52l) has no leading-twist contribution. This is also the reason, why the 
relative error in (|6Q|) is much larger than in (l59l) . The agreement of both results is still very 
reasonable, taking into account the smallness of Rem that is largely due to cancellations. It is 
therefore not unexpected that two different sum rules for the same quantity agree only within a 
factor of two and do not contradict the validity of our approach. 

In ll33l the form factors of the nucleon-A transition have been calculated as part of a very 
general examination of radiative decays of decuplet baryons into octet baryons. This calculation 
has been carried out in LCSR and is, in principle, very similar to our calculation. Apart from 
some technical differences, such as a different choice for the Lorentz basis, there is one point that 
has to be addressed. We explicitly included the electromagnetic background field in the quark 
propagator (PT4l) . As discussed in [|27l . working in this background field simplifies the treatment 
of the notorious contact terms. In particular, this procedure allows to include in a natural way 
the contributions from photon distribution amplitudes, that are known to vanish exactly, but 
have a non-zero conformal expansion to next-to-leading order. These have to be taken into 
account as most photon DAs are only known to next-to-leading order in conformal spin, which 
requires all photon DAs to be of this accuracy. In ll33l such terms were neglected. The numerical 
impact of these contributions apparently is small and the final results for G*m(0) = 2.5 ± 1.3 
and R E m(Q) = —6.8% from Ref. ll33Tl are close to ours. This consistency lends support to the 
general technique of LCSRs using photon DAs. 



4 Conclusions 

We have calculated the nucleon magnetic moments and the electromagnetic transition form fac- 
tors of 7p — > A + for real photons using the light-cone sum rule approach. Our result for the 
magnetic moment of the proton fi P is in good agreement with experiment and lends support to 
the current models of photon DAs and the estimates of the magnetic susceptibility of the quark 
condensate. The sum rule for the neutron magnetic moment is more sensitive to higher twist 
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terms and therefore less accurate. 

The calculation of the 'jp — ► A + transition form factors is, in principle, analogous to the 
calculation of the nucleon magnetic moments. The main difference is the asymmetric choice of 
the Borel parameters, that allowed us to take the mass difference between proton and A + into 
account. Refraining from doing so would lead to distinctly worse results. Within uncertainties, 
the magnetic dipole form factor Gm(0) of the jry — > A + transition agrees well with current data. 
This result is rather surprising as a different approach 113511 also using light-cone sum rules, that 
is valid for Q 2 > 1 GeV 2 predicts a value for Gm that is below data in the region Q 2 < 2 GeV 2 . 
In order to close the gap to this calculation it is necessary to expand our approach from the 
real photon point to virtualities ranging from to — lGeV 2 . This requires photon distribution 
amplitudes for virtual photons, see e.g. Il53l . 

Our results for the ratio Rem agree with experiment within a factor of 2. The two different 
sum rules written down for Rem also differ from each other by a factor of 2, while they agree 
very well for G*m(0). This supports our presumption that a lower accuracy for Rem is due 
to considerable cancellations, so this quantity is intrinsically more difficult to calculate with 
precision. Both G M (0) and Rem(Q) are in good agreement with the corresponding results from 

The best agreement for fip with experiment could be archived by the choice 

<p(l/2) = 5.25±0.15GeV~ 2 

and this value would also be favoured by the two sum rules for Gm(0), provided the asymptotic 
shape for (p(u) is used. The value of (p(l/2) is an independent piece of information compared to 
the expansion in Gegenbauer polynomials, which can be obtained from the sum rules. However, 
this is not sufficient to provide insights into the shape of the distribution amplitude. 

In order to increase the accuracy of our calculation it is necessary to take into account a s 
corrections. As the main source for the uncertainties are the numerical values of the twist- 3 and 
twist-4 DAs that are known only up to at best 50%, it is an important task of its own to determine 
their values more precisely. This would improve the numerics in this paper and be valuable for 
future work. 



Acknowledgements 

I am very grateful to V. Braun for many enlightening and clarifying discussions and to A. Lenz 
for proofreading the manuscript. I also would like to thank T. Aliev for pointing out the work of 
Ref.ES. This work is supported by DFG - 9209070. 



19 



A Photon distribution amplitudes 



For completeness we collect the relevant photon distribution amplitudes for the pj — > A tran- 
sition according to [|27l . Note that in 11271 the photon momentum has the opposite sign and the 
parametrisation of the separation of antiquark and quark is different. 
The path-ordered exponents (cp. Eq.(fT5l)) 

[x, y] = Pexp |z Jdt(x- y)^ [e q A»(tx - ty) + gB"(tx - ty)] 

assure gauge invariance of the matrix elements. It is important, that the electromagnetic field is 
included herein, as additional terms to those given in [27 J will occur otherwise. 



A.l Twist-2 and Twist-4 DAs 

The leading-twist DA reads 

(0| q(0) [0, x) a aj3 q(x) |0) F = e q (qq) I du ip(u)F a/3 (ux 



with 



+ 



16 



du x 2 A(u)F a p(ux) 



g (qq) r 

H — ^— — / du H>(u)x p (x/3F ap (ux) — x a F/3 p (ux)) 
° Jo 



ip(u) = ^ asy -(w) = 6xu(l- u) 

A(u) = 40u(l - u) (3k - k + + l) + 8 (C 2 + - 3C 2 ) x 

[u(l -u)(2+ 13u(l - u)) + 2u 3 (10 - 15u + 6m 2 ) ln(w) 
+2(1 - u) 3 (10 - 15(1 - it) + 61 - u 2 ) ln(l - u) 



B(u) = 40 f da(u- a) (l + 3k + ) 
Jo 



(62) 



(63) 



(64) 
(65) 



(0\q(0)e q [0,x] F fiv (ux)q(x) \0) F = e q (qq) / Pa<S 7 (a)F Mi/ (a u s 

(0|g(0)e g {0,x]a a pF^(ux)q(x) \0) F 
e q (qq) 



qx 



[qaq^xp - qpq^e^ - q^e^x p + g^e^zj 7?(u, qx) 



(0\q(0) [0,ux]gG IM/ (ux) [ux,x]q(x) \0) F = e q (qq) VaS(a)F )lu (a u x) 



(0|g(0) [0, ux\ i'fsgG^ (ux) [ux, x] q(x) \0) F = e q (qq) / VaS(a)F /J , 1/ (a u x 



(66) 

(67) 
(68) 
(69) 
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(0\q(0) [0, ux\ (T af3 gG^ u (ux) [ux, x] q(x) |0) ; 



Here we used 



and 



with 



-e, {qq) 


(A) 1 (A) j 


(A) I (X) 1 

q<*exu9^ + qf}e±v9aii 


Ti(u,qx) 




- e q {qq) 


(A) _L (A) _L 


(A) J_ , (A) X 
Qv e ±a9/3fi + qv e ±/39afi 


T 2 (u,qx) 




e q {qq) 
qx 


[q*q»e%xv - qpq»e± a - 


■ qaque^px^ + q$q v e\ a 


x^] T 3 {u,qx) 




e q {qq) 
qx 


[qaq^e^xp - qpq^e\ v - 




x a \ Ji(u, qx) 


(70) 


p 

/ Va 
■j 


= / d<y q / da-q 
Jo Jo 


I da g 5(1 — a q — a-q 
Jo 


-OCg) 


(71) 




= a q + uotg 






(72) 


9/j,i/ 


q^v qu^u 

= 9v* 

qx 




(73) 




= <7,V (A) 






(74) 



5(a) = 30a 2 [(k + «+) (1 - + (& + d + ) (1 - a s ) (1 - 2a 9 ) 

+C 2 (3 {a-q - a 9 ) 2 -a g (l- a g ))] 
5(a) = -30a 2 [(« _ «+) (1 _ Qp ) + ( Cl _ £+) (1 - (1 _ 2a 9 ) 

+C2 (3 (a f - a 9 ) 2 - a ff (1 - a g ))] 
5 7 (a) = 60a 2 , (a, + a g ) (4-7 (a 9 + a,)) 



7^(M,gx) = y Dae !a » ?1 T,(a) 
?i (a) = -120 (3C 2 + C 2 + ) («? - "g) ttgagttg 

r 2 (a) = 30a 2 {aq - a q ) [(« - k + ) + (Ci - Ci + ) (1 - 2a 9 ) + C2 (3 - 4a 3 )] 



Q a g 



T 3 {a) = -120 (3C 2 - C 2 + ) (a? - ««) «g« 
T 4 (a) = 30a 2 (a ? - a 9 ) [(« + « + ) + (Ci + Ci + ) (1 - 2a 9 ) + C2 (3 - 4a 9 )] 
Tj{a) = 60a 2 (a, - a q ) (4-7 (a 9 + at,)) . 



(75) 

(76) 
(77) 

(78) 

(79) 
(80) 
(81) 
(82) 
(83) 



The abbreviation a represents (a q , a-q, a 9 ).The values of the various constants can be found in 
tabled! 
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It should be noted that the matrix element 

(0| q(0)e q [0, x] cr a pF fJ/l/ (ux)q(x) |0) F 

vanishes exactly if one sums up the whole conformal expansion. The expansion itself has, 
however, non-zero coefficients and thus in next-to-leading order in conformal spin the matrix 
element is different from zero. For the same reason the matrix element 

(0|g(0)e, [0,x}F^(ux)q(x) \0) F 

has herein mentioned form and not e q {qq) F pv (ux) . 

A.2 Twist-3 DAs 



Where 



(0|g(0) [0,x]>y a q(x) \0) F = -| fa J du \u)x»F pa 
(0| q(0) [0,x] lal5 q(x) \0) F = -i jfa J du^ {A \u)x p F pa 



(0|g(0) [0,ux] igjaG^ux) [ux,x\q(x) |0) 7 



e qfaqa 



q c ±fl q e ±l/ 



q(0) [0,ux] g^alsG^iux) [ux,x] q(x) \0) F 

— e qfaqa 



q e ±n q e ±u 



^ (V \u) = -20u(l - u)(2u- 1) 
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+ — K " 3<) u(l - u)(2u - 1) (7(2u - l) 2 - 3) 



^ A \u) = (1 - (2u - l) 2 ) (5 (2u - l) 2 - 1) | (l + 



9 v 3 A 
16^ - 16^ 



V(a) = 540c^ (a q — a q ) a q a q a 2 



A(a) = 3Q0a q a q a g 



1 + ^(7^-3) 



(84) 
(85) 

(86) 

(87) 



(88) 

(89) 
(90) 
(91) 
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.3 Numerical values for the 

lGeV 



parameters at the renormalisation scale ji = 



X 


3.15±0.3GeV" 2 


K 


0.2 


K + 





Ci 


0.4 


ct 





C2 


0.3 


ct 







-(4 ±2) ■ 10- 3 GeV 2 




-2.1 ± 1.0 




3.8 ±1.8 


(9?) 


-(240±10MeV) 3 



Table 1 : Numerical values and uncertainties of the relevant parameters E71 l25l . 



Light-cone expansion of the —> A + correlation function 

jre we collect the results for the Feynman diagrams in Figj2]up to twist-4 accuracy. 
• Diagram a: 



T a (q,p) 



Diagram b: 

Tb(q,p) ~- 



7T- 



-up\ — up\ 



H (up 2 + 3u 2 p 2 + 8u 3 qp) ln ( 

64 v ' \—u 



up 2 — up 2 



(3p 2 + 3up 2 + 6uqp — 6u 2 qp) ln f — — 

64 v ' \—u 



up\ — up\ 



64 



uu 



In 



-up 2 — up\ 



■n, 



H (3p 2 + up 2 + Auqp + Au 2 qp) ln f — = 

64 ' \—u 



up\ — up 2 



■Tic 



+— [up 2 + 3u 2 p 2 + 8u 2 qp) ln 



-up I — up 1 



■TZ 



12 



+ 



' d (qq) / du 



(92) 



Up{ — UP2 



23 



2 I - K 3 



1 ( u )l n ( ^ ) A ^ 1 



Id — — 1^2, 



+ (^ ( _, p? _^ )ln (_^_) 



16 y-upi-up^ 



24 y — up\ — 



^ A \u) + 2i> { y\ u ) i 

~ u 75 — 2 ' n & 

48 — -up{ — -up2 



-^ A \u) + 2^ v \u) u f /i 2 



+ ~ v —In 



48 \ —up\ — up\ 



^ A \u) + 2^ v \u) 

24 —up[ — up\ 



-uqp — — 2 ) • K 9 



i boo 1 

H — 77; =~1 2 ' 

lb — up{ — up^ 

V 6 \-up\-up Z 2j 8 -up\-up z 2 J 
u^ A \u) + 2u^ v \u) ^ ( /i 2 



+ v y : r ^ln 



8 y — up\ — up\ 

u 2 qp{^ A \u) + 2^ v \u)) 1 



=~1 ^2 I ' ^12 



12 — up} — 

+ ... . 



(93) 
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• Diagram c: 







x { ( -]Sy(a)—-^ -r - ^—^V(a)—-^ -A ■ 1Z±+ 

1 4 J -(p + a u q) 2 4 -(p + a u q) 2 J 

+ f( S7 ( 2 ) + 77( S ))l„(- ¥ ^).7 J6+ 
■u 1 

+ 2^te) + gte)) _ (p+n>9) , -R M 

+ . . . . (94) 

• Diagram d: 

T d (p, q) = (p, q) + Tg ] (p, q) + (p, q) + T< 4) (p, q) + T d (5) (p, q) + . . . . (95) 
With 

>(p, *) = (e u - e d ) (qq) j\u j 'j>a S(a)±{ - • ^ 

+ (1-2m) " u -ft 5 
-Q» + a u g) 2 

(96) 



87r 2 -(p + a u q) 2 4 
(1 - 2w)a u 1 

• /C5 



87r 2 -(jo + a u g) 2 
Svr 2 n {-(p + au q) 2 ) 4vr 2 -(p + a n g) 2 
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(1 — 2u)a u 1 

2tt 2 ~(p + a u qy 



Tin } (97) 



T^(p, q) = (e u - e d )/ 37 J du \Va V(a) 



u 



qp 



+ 



In 



+ 



Air 2 -(p + a u q) 2 

u 1 

87r 2 -(p + a u q) 2 

ua u qp 



+ 



■K 3 



An 2 \-(p + a u q) 2 J 2ir 2 -(p + a u q) 2 



■11 



12 



(98) 



T?\p, q) = (e u - e d )f 3l J du j Vn A(n ) 



u 



\ An 2 


-(p + a u q) 2 


u 


1 


8tt 2 


-(p + a u q) 2 


u 


qp 



+ 



In 



An 2 \-(p + a u qf 



+ 



2n 2 -(p + a u q) 2 

(1 - 3u)a u qp 
2tt 2 -{p + a u qf 



■K 3 

■Kg 

n l2 



(99) 



Tf\p,q) = -(e u -e d ) (qq) [ 

Jo 

x <! / Va 



du 



u 1 — 2u u \ 1 

87T 2 + ~ 8^ T ^ ) -{p + a u q) 2 ' Ka 



+ 



+ / ©a 



ua u qp ua u qp , 
22(a) r-^TJa, 



+ 



u 



Air 4 



Air 4 



2tt 2 
1 



(>) _ 7(4) j 



+ 



2tt 2 



1 5 — 6m 1 — 2-u 3 — 2u 



8tt 2 



4tt 2 



x In 



8tt 2 
/i 2 



-(p + a u qf 



■Tift 



+ 



u 

87 2 



u 



(«)) 



7e 



10 
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+ 



where 



Va 



ua Ur ^ , . (1 — 2u)a Ur ^ , . wa Mm . . 
^T 2 (a) + 1 -V " T 3 (a) - -^T 4 (a.; 



2tt 2 



+ - 



2tt 2 



J(D + J(4)A + 



7T- 



2tt 2 
1 - 2m 



-(p+ a u q) 2 



71 s 



j(2) + j(3) 



ii 



Pa = / da q / dag / da g 5(1 — a q 
Jo Jo Jo 



a u = a q + ua g 



P 1 ' := u I da,, I da' 



l-a' 



da 



— u da' 



l-a!„ 



+ u da' 



da-n 



del' 



l-a'J, 



q ~(P + ( u + ua q - ua q )q) 2 

, Tija'^a^l- a' q - a' q ) 
« -(p+(l-n^)) 2 

lifai', a' 1-a'l-aL) 



da- 



da'-^ 



■(p + (u + ua')) 2 



(100) 

(101) 
(102) 



(103) 



Here fj, 2 is an arbitrary scale and the dots denote polynomials in p 2 and p 2 , that will vanish 
after a Borel transformation. The /W arise due to partial integration of terms proportional to --k 

The functions (p(u), A(u), B(u), 5(a), 5(a), <S 7 (a), V (A) (V)> V(a), .4(a), T^a) and 

7^ 7 (a) are defined in Appendix lAl 



C The functions X\, J 2 and X 3 



In this Section we have gathered the explicit expressions for the three functions I l5 X 2 and X 3 
that appear in Eqs.(l50l)-(l52l). 



2i=- 



t>A(v 



) f ^ _ e -So/t) 



+ 



da„ / da~a~^~^ — — <S 7 (a,j, a^, 1 — a q — a q ) 



o Jo 



2(1 — a q — a q ) 2 



v rv 



a I a I -%v + a q - a q ~ 
da n I da„ v—rrz — £> (a 9 , a q} 1 — a q — a q ) 



o Jo 



2(1 - a, - a ? ) 2 



ii rv 



+ / da / da-g v 



l-2v + a q - a q 

-1% a„, oc-n, 1 - a„ - as 



</0 



./0 



da„ / da-n 



2(1 — a 9 — a q 

v 



— J-3{W q , U> q , L — Ui q — Ut q ) 



2(1 - a q - a q 



■T 4 (a g , a ? , 1 — a q - a q ) 
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'0 JO 



— / da„ I da-aV 



+ 



I da q / da' q I da 
Jo Jo Jo 



2(1 — a q — OLq) 

v — a 



[T 2 - T 4 ] (a q , a q , 1 - a q - a q ) 



(l-a q - a- q 



x [-Ti + T 2 + T 3 - T 4 ](a' q , a q , 1 - a' q - a- q 



+ / da a / da' r 



l-a' 



da- 



v — a n 



q {l-a q - a q ) 2 
x [— Ti + T 2 + T 3 - T 4 ](a' q , a q , 1 - a' q - a q ) 



pv r l - a ' q 
\ da> q 

JO Jv 



da-H / da^r 



, v 



q a 2 



x [-Ti + T 2 + T 3 - T 4 ](a' q , a' q , 1 - a' q - a' q ] 



v pa 



+ da' da'' 



o Jo 



l-a' n 



da'-, 



q (1 - 

x [-7\ + T 2 + T 3 - T 4 ](a q , aL, 1 - a' q - ai) 



2 / da g / da' / cfag — — r%- 

Jo Jo Jo [l-a q -a q y 



2 / da„ / da' 



l-a 



x [T 2 + T 3 ](a;,ag,l-a;-ag) 

9 da- — - ^ — - a ^ 
" q (l-a q - a q Y 

X [T 2 + T 3 ](«g, ttg, 1 - a^ - ttg) 



+ 2 / da' 



a « , v 2 



da q I daL ^ [T 2 + T 3 ] (a/ , ai, 1 - a£ - a^) 



o 



-2 da' q da" / dai — [T 2 + T 3 ](aJ, 1 - aj - aQ 

JO JO JO U ~~ ^gj 



+ / da c 



o JO 



1— u 



da-n v 



1 — v + a„ — a 



-T2(a q , a q , 1 - a q - a q ) 



q 2(1 - a q - a q ) : 



x t (1 - e~ 5o/ *) 



(104) 



T9 = - 



/s 7 f v 



8tt 2 V3 

/3 7 



87T 2 7 



u r v - a 

da n I da, 



9 (I - a q - a q ) 2 



x [V -A] (a q , a q , 1 - a q - a q )t (l - e' Si)/t ) 



(105) 
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da„ / da„ v 



v — a„ 



q " (1 - a q - CXq) 2 
x [<S 7 + T 4 7 ] (otg, a q , 1 - a, - of J )t(l - e" 5o/ ') 



1 — 2d + a a — a 5 

9 °' 



- / da„ I da^v ' — - — 
2(1 - a q - a q ) 2 



x 



S + S 



(a q , acq, 1 - a q - a q )t(l - e 5o/ *) 



— / da q / da^D 
Jo Jo 



v — a T . 



2(1 - a q - a q ) 



x [T 2 - T 4 ] (a ? , a q , 1 - a q - a q )t(l - e^') 



+ 2 / da / da' / da 



(v - a Tl 



(1 - a„ - a 5 



x [T 2 - T 4 ](«;, aj, 1 - al - a,)i(l - e- So/t ) 



+ 2 1 da a \ da' 



l-a'„ 



da 



q (1 - Qfg - ttg) 3 



x [T 2 - T 4 ](a' a,, 1 - q/ - a f )t(l - e^') 



-2 / da' 



l-a' 



das 



? da;" K r^ 



<y-r, 



+ 2 / dal 



x [T 2 - T 4 ](a;, 4, 1 - a' q - a£t(l - e^) 



da„ 



da^ 



(l-< 



x [T 2 - T 4 ](a" a' 1 - a'' - a' )t(l - e"*/*) 



(106) 



The functions y>(u), A(u), B(u), 5(a), 5(a), <S 7 (a), ^ A )(w), ^ (y) (V>> V(a), .4(a), T;(a) and 
X±(a) are again defined in Appendix lAl 
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